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Abstract 

The impurities of exchange couplings, external magnetic fields and Dzyaloshinskii-Moriya (DM) inter- 
action considered as Gaussian distribution, the entanglement in one-dimensional random XY spin systems 
is investigated by the method of solving the different spin-spin correlation functions and the average mag- 
netization per spin. The entanglement dynamics at central locations of ferromagnetic and antiferromagnetic 
chains have been studied by varying the three impurities and the strength of DM interaction, (i) For ferro- 
magnetic spin chain, the weak DM interaction can improve the amount of entanglement to a large value, 
and the impurities have the opposite effect on the entanglement below and above critical DM interaction, 
(ii) For antiferromagnetic spin chain, DM interaction can enhance the entanglement to a steady value. Our 
results imply that DM interaction strength, the impurity and exchange couplings (or magnetic field) play 
competing roles in enhancing quantum entanglement. 
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I. INTRODUCTION 



Entanglement is not only important in the quantum information processing (QIP), such as quan- 
tum teleportation,' 1 ' dense coding, ' 2 ' quantum secret sharing, ' 3 ' quantum computation,' 4 ] but also 
relevant to quantum phase transitions' 5 ' in condensed matter physics. In order to realize quantum 
information process, great effort has been devoted to studying and characterizing the entanglement 
in cavity QED.' 6 8 ' Now, much attention has been paid to the entanglement in spin systems, such 
as the Ising model ' 9 ' and all the kinds of Heisenberg XY XXZ XYZ models. ' 10 ~ 13 ' However; as 
far as we know, most discussions mentioned above merely focused on the models with spin ex- 
change couplings, while Dzyaloshinskii-Moriya interaction has seldom been taken into account. 
The antisymmetric DM interaction, introduced by Dzyaloshinskii and Moriya, is a combination of 
superexchange and spin-orbital interactions. In fact, some one-dimensional and two-dimensional 
spin models have manifested such interactions J 1415 ] Therefore, it is worthwhile including DM in- 
teraction in the studies of spin chain entanglement. 

Impurities necessarily exist in real materials and their effects are more pronounced in con- 
densed matter physics. Thus, it is important to study the effects of impurities in view of the pos- 
sible realizations of one-dimensional ferromagnetic and antiferromagnetic chains. In the previous 
researches, the impurity effects on the quantum entanglement have been studied in a three- spin 
system' 16 ' 17 ^ and a large spin systems under zero temperature.' 18 ' However, in these works, they 
have just studied single impurity. 

Recently, Huang et a/.,' 19 ' Osenda et a/.' 20 ' and we' 21 ' have demonstrated that for a class of 
one-dimensional magnetic systems, entanglement can be controlled and tuned by introducing im- 
purities into the systems. For the pure case, Osterloh et al. ' 22 ' examined the entanglement between 
two spins of position % and j in the spin chain as the system goes through quantum phase transi- 
tion. They demonstrated that entanglement shows scaling behaviour in the vicinity of the transition 
point. For a two-qubit spin chain with DM interaction, researchers' 23 ' 24 ' have considered thermal 
entanglement and teleportation. For a particular spin system the allowed components of the DM 
interaction are determined by the corrections to the energy symmetry of the spin complex. Since 
the DM terms break spin- spin rotational symmetry, we need to calculate how spin exchange cou- 
plings and DM interaction have effect on the entanglement and phase transition point. It is an 
interesting quantum phenomenon that the entanglement shares many features with quantum phase 
transition (QPT), QPT is a critical change in the properties of the ground state of a many body sys- 
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tern due to modifications in the interactions among its constituents. The associated level crossings 
lead to the presence of non-analyticities in the energy spectrum. Therefore, the knowledge about 
the entanglement, the non-local correlation in quantum systems, is considered as the key to under- 
stand QPT. That is the purpose and motivation of the present work to investigate the behaviour of 
entanglement at and around the quantum critical point in one-dimensional XY spin system with 
DM interaction, which can display a variety of interesting physical phenomena providing new in- 
sight in two-site entanglement and the related QPT as well under the effect of the impurities of 
exchange couplings, external magnetic fields and DM interaction. 

We consider Heisenberg XY model of N spin-| particles with nearest-neighbour interactions. 
In the presence of impurities and DM interaction, l 25 ' one-dimensional Hamiltonian is given by! 19 ' 



i , N 1 N 
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where J^i+i and D iji+ i are exchange interaction and DM interaction along ^-direction between 
sites % and % + 1 respectively, hi is the strength of external magnetic field on site i, a x ' y ' z are the 
Pauli matrices, 7 is the degree of anisotropy and TV is the number of sites. For all the interval 
< 7 < 1 and N = 00, they undergo a quantum phase transition at the critical value A c = 1. 
The periodic boundary conditions satisfy <?% +1 = of, u y N+1 = erf, o z N+1 = a\. Let us define the 
raising and lowing operators af, a~ and introduce Fermi operators c+ and Cj,' 26 ' the Hamiltonian 
has the form 

N 

H = - Y^[(( J i,i+i + iD iji+1 )c+c i+1 + h.c) + 
i=i 

N 1 

(j i>i+11 44 +1 + h.c)] - 2 H4a - 3) ( 2 ) 

i=i 

In this study, the exchange interaction has the form Jj^+i = J(l + a iji+ i), where a introduces the 

impurity in a Gaussian form centered at ^ ^ ~ w i m strength or height (, a iji+ i = ( exp ^ — e (i — 

N + l\\ 

e is the value of the width of the distribution. For J < 0, the spin chain is antiferromag- 
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netic; for J > 0, the spin chain is ferromagnetic. The external magnetic field and Dzyaloshinskii- 
Moriya interaction take the form hi = h(l + fa) and D i>i+1 = D(l + 77^+1) exp (7^)' where 

fa = £exp [-e(i- — y— ))» V^+i = Kex P ( ~ ~ When « = = ^ = 0, 



we recover the pure case; when rj = 0, we recover the case described in Ref. [19]. For the dis- 
tributions of exchange interaction impurity, Dzyaloshinskii-Moriya interaction impurity and the 
magnetic field impurity, we fix the value of width of the distribution at e = 0.1 in all the cal- 
culations. As the center ^ ~ ) an( * tne w ^ m ( e ) °f me Gaussian distribution are fixed, we 
can obtain different impurities a, (3, r\ of the Gaussian distributions only by changing strengths or 
heights C, £, k. By introducing the dimensionless parameter A = J/2h, the symmetrical matrix A 



and the antisymmetrical B, the Hamiltonian becomes 



N 1 

H =J2 tfAijCj + M ih.rj + h.c) 



(3) 



i,j=l 



The above Hamiltonian can be diagonalized by making linear transformation of the fermionic 
operators r] k = 9kiCi + h ki cf, rj^ = J2? =1 9ki4 + h ki c u then the Hamiltonian becomes 
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H = ^2 A kVkVk + const, 



(4) 



k=l 



two coupled matrix equations satisfy 4> k (A — B) = A k ip k , 4>k(A + B) = A k (p k , where the compo- 
nents of the two column vectors 0^, ipki are given by <f) ki = Qh + h ki , ip ki = g ki — h ki . Finally, the 
ground state of the system \ip ) can be written as rj k \^ \ = 0. 

Using Wick's theorem, l 27 ' spin- spin correlation functions for the ground state and the average 



magnetization per spin can be expressed as 
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,. , — — J2k 'tpki'fikj Next, we give the expression of concurrence that quantifies the amount 
of entanglement between two qubits. 

For a system described by the density matrix p, the concurrence C reads^ 28 ! 



C(p) = max(0, Ai — A2 — A3 — A4) 



(5) 



where Ai, A 2 , A 3 , and A 4 are the eigenvalues (with Ai being the largest one) of the spin-flipped 
density operator R, which is defined by R = y^/pp^/p, where p — (a y <g> a y )p*(a y ® a y ); 



p denotes the complex conjugate of p; a y is the usual Pauli matrix. Using the operator ex- 
pansion for the density matrix and the symmetries of the HamiltonianJ 29 ' in the basis states 
{| TT)i I t!)> I It); I II)}' P has the general form. We can express all the matrix elements in the 
density matrix in terms of different spin- spin correlation functions. 

In this study, we focus our discussion on the transverse Ising model with 7 = 1. We examine 
the dynamics of entanglement in varying the impurities of exchange couplings, external magnetic 
fields and Dzyaloshinskii-Moriya interaction. First, we examine the change of the entanglement 
for the nearest-neighbouring concurrence C(i,i + 1) for different values of the impurity as the 
parameter A varies. Figure 1 depicts the nearest-neighbouring concurrence C(49, 50) as a function 
of the reduced coupling constant A at different values of the exchange couplings impurity ( and 
external magnetic fields impurity £ with the system size N = 99. Figure 1(a) shows the change 
of concurrence C(49, 50) as a function of A for different values of exchange couplings impurity 
with D — 0, i.e. in the absence of DM interaction. For the case of A > 0, we can see that the 
concurrence increases and arrives at a maximum close to the critical point A c , while it is close to 
zero above A c . As ( increases the concurrence tends to increase faster, and A m , where concurrence 
approaches a maximum, shifts to left very rapidly. This is consistent with the result reported in 
Refs. [19,21] (Fig.l). A similar behaviour can be seen for the case of A < 0, that is to say, the 
entanglement has equal value for ferromagnetic and antiferromagnetic chains with the same | A | . 
The effect of the external magnetic field £ in the Gaussian distribution is also shown in Fig. 1(b). 
However, different from the effect of the exchange couplings, the concurrence increases slowly 
and tends to moving to infinity by increasing the value of the parameter £. This is also consistent 
with the result in Refs. [19,21] (Fig.l). In Figs. 1(c) and 1(d), taking DM interaction into account, 
we give a plot of the concurrence against exchange couplings impurity and external magnetic field 
impurity with D = 0.5 \J\. As ( increases, the concurrence increases slowly and the peak value 
decreases, which can be seen in Fig. 1(c), different from the result in Fig. 1(a) for ferromagnetic 
spin chain. Moreover, some interesting physical phenomena occur for the antiferromagnetic 
chain, for example, the concurrence decreases to zero at the critical point (A ) and increases from 
zero to a finite steady value across the transition point. Therefore, we can further understand the 
relation between the entanglement and quantum transition. In Fig. 1(d), the numerical calculations 
show that the steady concurrence decreases with the increase of £ for the antiferromagnetic 
chain, which indicates that the behaviour is very different from those in Fig. 1(b). Now weexplain 
why the curves of concurrence have some maximum or minimum at some special values of the 
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DM interactions and external magnetic fields. As [dC(4Q, 50)/d\]\ c diverges, the maximal or 
minimal entanglement will not occur at the critical point but in the vicinity of the transition point 

([<9C(49,50)/<9A] Am = 0, [<9 2 C(49,50)/<9 2 A] Am < or [<9 2 C(49, 50)/<9 2 A] Am > 0). In our model, 
when D = 0, quantum transition point A c = - ; from the above expression, we can know 

1 + <Xi,i+l 

the transition point shifts and is affected by the impurities of exchange couplings and external 
magnetic fields. For the entanglement length (or the correlation length), the position of the related 
maximal and minimal concurrence will shift in the same way. However, DM interactions lead 
to different coefficients in the first two parts of Eq.(2), so the critical point A c occurs between 

the two ones (\^- — ^±1, — -V For the case of J < 0, there exist two 

\l + a iji+ i 1 + a i:i+1 - + rji ti+1 )/ 

critical points. Of course, we can figure out exact critical value and maximum or minimum of 
the concurrence through solving the first order and second order derivative of the entanglement 
respectively. 

From Fig. 1 , we can see that DM interaction plays an important role in enhancing entanglement, 
so it is necessary to study the effect of DM interaction on the entanglement. In Fig. 2, we show 
the results of the nearest-neighbouring concurrence as a function of the parameter A for DM 
interaction impurity at different strengths of DM interaction D. We can easily find that the 
competing roles played by DM interaction impurity k, strength D and exchange couplings J (the 
external magnetic field is fixed) in enhancing quantum entanglement will exist in spin chain. The 
competing effect leads to shift of the critical point and the entanglement. The results show that 
when the absolute value of A is below A c , the concurrence only increases with |A|, DM interaction 
impurities will have no effect on the entanglement once strength is fixed, i.e. exchange couplings 
is predominant in the competing role. The effect of weak DM interaction strength D = 0.1 |J| 
is shown in Fig. 2(a). Contrast to the exchange couplings impurity, DM interaction impurity 
can enhance the entanglement, the concurrence increases and tends to move to infinity(> 0) by 
increasing the value of the parameter k. It is interesting to find that the entanglement peak and 
steady value between the nearest neighbours with D = 0.5 \J\ increase to a value larger than 
those in Fig. 2(a). With the increasing D, in Fig. 2(c), the concurrence decreases as k increases. 
We can imagine that there must be a critical DM strength (D c ), below D c , impurity enhances 
entanglement, while above D c , impurity shrinks entanglement. In other words, at some special 
values of the DM interactions, the entanglement varies at different critical vicinities, which is 
similar to the analysis in Fig.l. The comparison among the different curve in Fig. 2(d) shows that 
the concurrence decreases rapidly above A c by increasing the value of the parameter k, which is 
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different from the results obtained from Figs. 2(a) and 2(b). That is to say, the strong D is not 
helpful to keeping the better entanglement for Gaussian distribution. 

The effect of DM strength is demonstrated in Fig. 3 by the evolutions of the concurrence. 
Figure 3(a) corresponds to the case of k — 0, the nearest-neighbouring concurrence increases 
with the increase of D, a critical point occurs with small DM interaction strength and the peak 
of the maximal entanglement becomes larger. It is the DM interaction that leads to considerable 
different evolutions of the entanglement, hence the entanglement is rather sensitive to any small 
change with the DM interaction. Thus, by adjusting DM interaction one can obtain a strong 
entanglement. Similar behaviours to those in Figs.2(c)and 2(d) are shown in Figs. 3(c) and 3(d), 
we can see that DM interaction strength is not certain to enhance the entanglement, and the 
entanglement tends to be reduced in the presence of strong DM interaction at k = 1.0. The results 
we have obtained here are also consistent with those in Fig.2. 

According to finite-size scaling analysis, the two-site entanglement is considered as a function 
of the system size (including the thermodynamic limit) and the distance |A — A c | from the critical 
point. The entanglement can approximately collapse to a single curve for different system sizes 
ranging from 41 up to 401, thus all key ingredients of the finite-size scaling are present in the 
concurrence. The first order derivative around the critical point becomes sharper (the peak 
position A m in approaches the critical point A c ) as the system size increases, and is expected to be 
divergent in an infinite system (dC/dX = Ai In |A — A c | + const). Though there is no divergence 
when iV is finite, the anomalies are obvious. Its value diverges logarithmically with the increasing 
system size as dC/d\ = A 2 In N + const. Thus, we can see that the QPT of the system is 
reflected by the behaviour of the concurrence and its A derivative and finite size scaling is fulfilled 
over a very broad range of values of N, which are of interest in quantum information. 

In summary, from the above analysis, it is clearly noted that the three different impurities and 
DM interaction strength, which play the competing roles in enhancing quantum entanglement, 
have a notable influence on the nearest-neighbouring concurrence in the one-dimensional s = ^ 
random XY spin system. The nearest-neighbouring concurrence exhibits some interesting 
phenomena. For an antiferromagnetic spin chain, there is a critical point where the entanglement 
is zero. DM interaction is predominant in the competing role and can enhance the entanglement 
to a steady value. For a ferromagnetic spin chain, the weak DM interaction can improve the 
amount of entanglement to a large value. However, under condition of strong DM interaction, 
there is a critical point D c where the impurities have the opposite effect on the entanglement 
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below and above D c . Thus we can employ DM interaction strength as well as three different 
impurities to realize quantum entanglement control. For the case of 7 7^ 1 (XY model) or the next 
nearest-neighbouring concurrence related QPT, we will present further reports in the future. 
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FIG. 1: The nearest neighbour concurrence C(49, 50) for the impurities £ and the impurities £ as a function 
of the reduced coupling constant A with N = 99, 7 = 1, k = 0. 
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FIG. 2: The nearest neighbour concurrence C(49, 50) as a function of the parameter A for the DM interac- 
tion impurities k with N = 99, 7 = 1, £ = £ = 0. 
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FIG. 3: The nearest neighbour concurrence C(49, 50) as a function of the parameter A for the DM interac- 
tion strength D with AT = 99, 7 = 1, C = £ = 0. 
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